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Abstract. This paper studies the properties of the probabihty density func- 
tion Pa.f.nCx) of the n-variate generalized Linnik distribution whose charac- 
teristic function ipa.i/.n(t) is given by 

iPc,^,n{-t) = -, n'-n-TF. a e (0,2], 1/ > 0, t e R", 

(1 -I- lltll") 

where ||t|| is the Euchdean norm of t £ R". Integral representations of 
Pa,iy,n(x) are obtained and used to derive the asymptotic expansions of pci,i/,n(x) 
when ||x|| — ► and ||x|| — > oo respectively. It is shown that under certain con- 
ditions which are arithmetic in nature, Pa,i^,n(x) can be represented in terms 
of entire functions. 



1. Introduction 
In 1953, Linnik showed that the function 

(1-1) "^"(^^^TTl^' *eM, 

is the characteristic function of a symmetric probabihty density pa {x) . Since then, 
Pa{x) is known as Linnik's probabihty density. When a = 2, P2{x) is the probabihty 
density of the Laplace distribution. Therefore, Pa{x) is sometimes referred to as a- 
Laplace probabihty density. In |26^ , Klebanov, JVIaniya and Melamed introduced the 
concept of geometric strict stabihty and showed that the characteristic functions 
of the probabihty distributions that are geometricaUy strictly stable are up to a 
scaling factor given by: 

1 / TTQ; TTd \ 

(1-2) ^^At) = Y^—e7i^^ «G(0,2), |0|<min(— ,7r-— ). 

The corresponding probability distribution is known as asymmetric Linnik distribu- 
tion. In [64], Fakes showed that the probability distributions known as generalized 
Linnik distributions which have characteristic functions 
(1.3) 

1 / TTO. 7T(y.\ 

^^Mi)= (i^g-.,9sgnt|^|o)^ ' «e(0,2), 1^1 <min(^— j , v>0, 

play an important role in some characterization problems of mathematical statistics. 
The class of probability distributions with density functions Pa,v,9(x) have found 
some interesting properties and applications [2lltll5l[7ll8l[9llM | [29 t [30 l [32 l l33 l 
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[MlESlEelElllSHlllTlllalllllllSllS^ in particular, 

they are good candidates to model financial data which exhibits high kurtosis and 
heavy tails. In [69j and [T], the probability distribution Pa{x) was independently 
generalized to multivariate distributions. Anderson ^ showed that the function 

(1.4) (/?a,n.s(t) = ^, a e (0, 2], t e M", 

1 + [i'Et] 2 

where E is a positive definite n x n matrix, is a characteristic function of a n- 
variate probability distribution which he called the n-variate Linnik distribution. 
Properties and applications of this distribution have been studied in UHl ISZl 
l39l l40l |4T[ l42l l43l [60] . Similar to the univariate case, the multivariate Linnik 
distribution is a subclass of multivariate geometric stable distributions [SDl [321 ■ 

As the case of stable distributions, in general the probability density functions 
of the univariate and multivariate Linnik distributions do not have closed forms. In 
|28j . Kotz, Ostrovskii and Hayfavi studied the analytic and asymptotic behaviors 
of the probability density function of the univariate symmetric Linnik distribu- 
tion whose characteristic function is given by (jl.ip . This work was generalized to 
symmetric homogeneous multivariate Linnik distribution with characteristic func- 
tion (|1.4p with E = id by Ostrovskii [62j : and to asymmetric univariate Linnik 
distribution with characteristic function (jl.2p by Erdogan [13] , and finally to asym- 
metric univariate Linnik distribution with characteristic function (jl.Sp by Erdogan 
and Ostrovskii [Til [H] ■ In the present work, we consider the generalization of the 
multivariate Linnik distribution whose characteristic function is given by 

(1.5) iPc.,.,n.s{t) = —17, aG (0,2], z/>0, teM". 

(^1 + [t'Et]^j 

One can expect that this multivariate generalized Linnik distribution also plays an 
important role in some characterization problems in multivariate statistics as in the 
univariate case ^64J. The main goal of this paper is to investigate the asymptotic 
and analytic properties of the probability density function of the multivariate gen- 
eralized Linnik distribution. It suffices to restrict to the case where E = id, i.e., we 
only consider the class of multivariate symmetric generalized Linnik distribution 
with characteristic function given by 

(1.6) = (TTMPf ' aG (0,2], z^>0, « SK", 

where ||t|| = ■^/^"^j^ i| is the Euclidean norm of t. We denote by Pa_i/^„(x) the 
corresponding probability density function. 

There are also other motivations to the problem we want to study here. Recall 
that the dual distribution p of a probability distribution p is a probability distri- 
bution whose characteristic function (resp. probability density function) is up to a 
constant, the probability density function (resp. characteristic function) of p |15j . 
When a € (0, 2] and v > n/a, the dual distribution of the generalized Linnik distri- 
bution is called the generalized Cauchy distribution |25| whose probability density 
function is up to a constant, given by (|1.6p . Therefore, Pa,i^,ni^) /Pa,i^,niO) is the 
characteristic function of the generalized Cauchy distribution. For the special case 
where a = 2 and v > n/2, the generalized Cauchy distribution is the well-known 
Student's t distribution which has a lots of applications (see e.g. [3T] and references 
therein) . 
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From the perspective of stochastic processes, Gneiting and Schlater [18] intro- 
duced a new class of stationary stochastic processes called Gaussian field with 
generalized Cauchy covariance whose covariance function is given by (|1.6p . This 
stochastic model has two parameters a and v which can give separate characteriza- 
tions of the fractal dimension and Hurst effect. It has found applications in many 
modeling problems [19l|46l|47l|48l[58l[59l[67l[2Q]. In this context, the function un- 
der investigation Pet. iy_„(x) appears as the spectral density function of the stochastic 
model On the other hand, when a = 2 and ly > n/2, up to a constant, (|1.6p 
appears as the spectral density of the Whittlc-Matern random field [Ml [731 HM ■ 
This class of random fields has found wide applications especially in geostatistics 
[l5l[2Ql[22l[23ll5l[56l[5a[lIll22]-In[50],we generalized the Whittle-Matern ran- 
dom field to a random field whose spectral density is given by (|1.6p with a € (0, 2], 
V > nja and showed that it can provide a more flexible model for wind speed data. 
The function pQ.iy,„(x) we want to study in this paper then becomes the covariance 
function of the generalized Whittle-Matern random field. 

From the brief discussion above, one notes that the function appears 
in various places and plays different roles, stretching from the probability density 
function of generalized Linnik distribution and spectral density function of Gauss- 
ian field with generalized Cauchy covariance to the characteristic function of the 
generalized Cauchy distribution and the covariance function of generalized Whittle- 
Matern field. Thus, a detailed study of the analytic and asymptotic properties of 
Pa,y,n(x) is of strong interest and significant importance, particularly in view of its 
potential applications in physics, internet traffic and financial time series analysis 
and modeling. This is the task undertaken in the present work. 



2. Multivariate generalized Linnik distribution 

Gneiting and Schlather [18] asserted that using similar arguments as of [16] and 
references therein, one can show that the function ipa,u.n{'t) (|1.6p is a covariance 
function of a stationary Gaussian random field if and only if a G (0,2] and v > 
0. This implies that when a G (0,2] and v > 0, Pa,!^,n(x) > and therefore is 
also the probability density function of a distribution. In this section, we provide 
another argument to show that (pa.u,n{t) (|1.6|) is the characteristic function of a 
probability distribution. The following is a generalization of the result of Devroye 
[8] to multivariate case. 

Proposition 2.1. Given a G (0,2], > and E a positive definite matrix of rank 
n, let Sa,'n,T. be a symmetric multivariate stable random variable with characteristic 

function exp ^— [t'Tit] ^ ^ and let U^, be an independent univariate gamma random 
variable with probability density function 



u > 0. 



Then the characteristic function of the random vector l^a,v.n,Y. — Uv^a,n,Y, is given 
by 1131). 
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Proof. 



E{e'^-^) =^(^S j^e'C^°t S".".^ 

=E (exp(-;7,[t'St]^)) 

= J-^^°°.-exp{-.(l + [m]t)} 



du 



1 + [t'^ty 



□ 



We call the random vector Xq_^^„ whose characteristic function is given by (|1.6p 
an (ck,j/, n) Linnik random vector. As in the introduction, denote by pa_i,^„(x) its 
probability density function. In other words, Pa,i/,ra(x) is the unique function such 
that 

(2.1) 'fc.,^.n{t) = (^i^^^^o^Y = j^^ e*''-Va.^,«(x)d"x. 

Since ^Pa,v,n{t') is a radial function, Pq.i/,„(x) is also a radial function. Denoting by 
<la,v,n('i') the function 

Qa,v.n{f) = J'a,i^,n(x)|jj^jj_^ . 

Then Schoenberg's formula gives 

1 . 
^^■^^ (l + ||t||o)'- ^ (2^)^11^11^ J^(r||t||)g„,,,„(r)r^(ir, 

where J^{z) is the Bessel function of the first kind. In this paper, we want to study 
the properties of the function qa,u,n{T)- 

Proposition 12.11 gives a relation between the probability density functions of 
the generalized Linnik distribution and the symmetric stable distribution. More 
precisely, let S'q^„(x) be the probability density function of the symmetric stable 
distribution with characteristic function exp(— i.e. 

(2.3) e'"*"" = / e'*-''s„,„(x)d"x, 
or equivalently, 

(2.4) s„,„(x) = -i-/ e-*-e-ll*ll°d"i=M^ Tj (i||x||)e-*°ttcit. 

(27r)" 7r„ (27r) 2 2 

Then 

Proposition 2.2. Let a e (0, 2] and v > Q. The probability density function 
PQ.i/,n(x) has a representation 

I nv-n-l-u"- I ll^l 



w/iere Sa,„(r) = Sa,„(x)L 
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The properties of the function Sa,nir) have been studied in p^. In particular, 
it was shown that 

Sa,n{r) =0{r°) as r ^ 0, 

(2-6) . _ _ N 

Sa,n{r) [r " ") as r^oo. 

These show that (|2.5p is indeed well-defined for all a € (0, 2], > and r > 0. 

Proof. Making a change of variable, we can rewrite the right hand side of p.Sp as 

^ .---1 /Ilxl 



m Jo 



(2.7) — ^ / u'^-^-^e-"S„,„ ^ U^^ 



Taking the Fourier transform of (j2.7|) and using the definition of Sa,n{r) (|2.3p . we 
obtain 



1 

1 



gmit.x^,.-lg-«~^^^ (||x||) d"xdM 
^.-lg-«(l+||t|r)^y 



1 



(i + ii*ihr' 

This proves the assertion of Proposition 12.21 The interchange of the order of inte- 
grations can be justified using (|2.6p . 

□ 



As is observed in [T3], the identity (see e.g. pT] ) 

J,,(z) = ^J^(z) - J^+i(z) 

and (HiH) imply that 

It follows from (^3]) that 

(2.8) Qa,,.,n{'r) = -27rrg„,^,„+2(?'). 

Since qa,u,n+2{r) > 0, we deduce from (|2.8p that 

Proposition 2.3. q2.u,n{r), f G *s a monotonically decreasing function of r. 

Now we want to find other representations for qaM^nif)- A naive application of 
Fourier inversion formula to (|2.ip gives 

^ — ix.t 

(27r)« (1 + ||t|h)^ 



(2.9) 



X 2 



^^(l|x||t) „ 

t-dt. 



(27r)t 7o (1+i") 
However, since (see e.g. [2T, ) 

2 / TT/^ TT 



as z ^ CXI, 



is only valid when av > ^^^^ . In [JS] , we proved the following: 
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Proposition 2.4. For a € (0,2) and v > 0, qa,v,n{T) has the following integral 
representation 

(2 10) ga,i/,n(r) = „-2 „+2 I'm- / 7 ^u^du. 

2 2 TT 2 Jo f 1 + e^u°^ 



For a — 2 and v > 0, q2.u,n{f) is given explicitly by 

(2-11) q2,.Ar) = :T^+^£r4p^^f-.W- 

_ffere Kf^{z) is the modified Bessel function of the second kind. 
Since Kf^{z) is an infinitely difFerentiable function and 

(2.12) if,(z)=o(|zrl-llog^) as.-O, 
K^{z) ~0{e^^') as z ^ oo, 

(see e.g. [21]), Proposition 12.41 shows that qa,i^,n{r) is an infinitely differentiable 
function of r when r S (0, oo). 

When 1/ = 1, (|2.10p has been proved by Ostrovskii [5^. When n = 1, it was 
proved in [51I1H1 Ell [12] ■ Both f^J^i and ([^H]) can be proved by showing that 
they satisfy (|2.2p . However, to motivate these formulas, one can start from (|2.9p 
when ai^ > ^^y^- In this case, the formula (|2.9p can be evaluated explicitly when 
a = 2 which gives (|2.1ip . When a € (0, 2), we use the fact that (see e.g. [21]) 

i/«(z) = J^(z) + ziV^(z), 

and 

K,,{z)^-e^Hl'Hiz), 

where Hj^\z) is the Hankel's function of the first kind and N^(z) is the Bessel 
function of the second kind or called the Neumann function. A change of contour 
of integration from the positive real axis to the positive imaginary axis gives (j2.10p . 
For details, see [i5] . 

Using Proposition 12.41 we can prove a result stronger than the result of Propo- 
sition [2T3I When a = 2, using the formula (see e.g. [2T] ) 

(2.13) K,{z) = / / e-y^iy' - iT-'^dy, fi>--,z>0, 
we can rewrite (I2.1ip as 

q2 t^Jr)^ 1 / (y^ - l)'^-''e-''ydy 

when V < When = ^^y^, we have the explicit formula 

Therefore (—i)'^ qi'^l ni''^) > for all k ~ 0,1,2,... and q2,v,n[r) is a completely 
monotonic function of r G M+ when v < When a G (0,2), we rewrite (|2.10p 



2-7r^r(ii±i; 
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as 



(2.14) ga^^,„(r) 



If n = 1, this gives 



f 2 



2 2 TT 2 Jo 



1 + e— 



1 + e 2 



■Kn-2 [ru)u 2 c?M. 



(2.15) 



9a, 1^,1 (f) 



1 



oo sm 



1 + e 2 



1 + e 2 

When n > 2, we can use (|2.13p to transform (|2.14p to 
(2.16) 



Qa^iy^n (^) 



2»-%T^r(^) 
1 



sin l^j/arg ^1 + e' 2° 



1 + e~-u" 
y sin ^i/arg ^ 



1 + e 2 It" 







1 + e 2 



(y2 - 1) — e-^™(i2/dw 



(2/^ — U^) 2 \ rfy. 



Since sin (uarg (l + e^^u"^^ > if ai/ < 2, ((TTSl) and imply that qc,,v,n{r) 

is completely monotonic if a G (0, 2) and av e (0, 2]. The result is summarized in 
the following proposition. 

Proposition 2.5. The Junction qa^u, nil") is completely monotonic if a G (0,2) and 
av £ (0, 2] ora^2andiye (O, 2±i] . 

We remark that when a G (0, 2) and v = \, this result has been proved in [62] . 
When n = 1 , it has been shown in [TTl [12| . 

3. Asymptotic expansions of gQ,i/,„(r) 

In this section, we derive the asymptotic expansions of the function qa.v,n{f) 
when r ^ and r ^ 00 respectively. We discuss the cases a — 2 and a G (0, 2) 
separately. 

3.1. a — 2. When a — 2, the asymptotic behavior of K^{z) as z ^ 00 (see e.g. 
[U]) gives 

Proposition 3.1. When a = 2 and u > Q, the asymptotic expansion of the function 
qa,u,n(j) when r —^ 00 is given by 

q2,i^,n{r) 



E 



1 ^C-^+J--),.-^-: 



In particular, the large-r leading term of q2^,j.n{r) is given by 

92, !/,«(?') 



+0 (r^'-^-^e-'-) 



r 2 6 



2V+''7rVr(i/) 

For r ^ 0, the explicit series representation of K^{z) about z — Q (see e.g. [21]) 



gives: 
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Proposition 3.2. When a ~ 2 and v > 0, 

I. If J — v is not an integer, then the series expansion of qa^ii.yiif) about r — is 
given by 

_ 1 ^ ( (r/2)^J'+^"-" (r/2)^J 

'^2-'"^^) - 2".-^r(.) sin [n (I - .)] ^„ [jlT {j + . - jir {j + 

II. If ^ — V = I is an integer, then 



, ^ ('r/2ll'l-'+2-'' r /r\2 

+ (-1)'+^ E ^yiTTyr r (2) - ^(-^^ + ^) - '/^("^i + ^' + ^) 

Notice tliat tlie beliavior of q2.v,n{r) as r depends on wliether ^ — v \s 
an integer, li ^-^ — v in not an integer, q2.v,n{r) can be represented by the sum 
of two series, one is an absolutely convergent power series, and the other is the 
multiplication of r^'^"" with an absolutely convergent power series. When ^ — v 
is an integer, q2,u.n{i') can be represented by the sum of three terms, one is an 
absolutely convergent power series, one is the multiplication of 7-ii"n{2i/-rt,o} 
a polynomial, and the third one is the multiplication of Inr with an absolutely 
convergent power series. 

Corollary 3.3. When a = 2 and v > Q, the r ^ leading term of the function 
qa,u,n{f') depends on the value of v: 



n 



2 J J2v-n 



I.Ifv<-, then q2,^,n , 
II. Ifv > -, then q2,y, - ^ ^' 



r 



2' 2«7rtr(i/) 
///. Ifv = ^, then C72,.,„ ^^J^ |- log (0% 2^(1)| . 

3.2. a e (0,2). When a e (0,2), the asymptotic expansion of qa^v.nir) as r — + 00 
can be obtained from the representation of qa,v.n{i') given by (|2.10p . 

Proposition 3.4. When a £ (0, 2) and ly > 0, the following formula is valid: For 
r 00, 



7^w(r), 



where 



2-^ T{N + u)T{^)T{^ ) 
n'^NWii^) sin''+^^ 



(3-2) \nNir)\ < ^ V 

2 
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In particular, the large-^ leading term of qa,u.n{r) is given by 



Q_a^v.n (^) 



■ sin — r 



n + 2 — r, 
TT 2 ^ 



r 



Proof. Making a change of variable in (|2.10p . we have 



(3.3) qc,i,Ar)=- 



n + 2 
2 2 TT 2 



Im 



{l + e—r-°'u<' 



^2 du. 



Using integration by parts, one can check easily that for an infinitely differentiable 
function g{y), we have 



{y-vr-'g^''Hv)dv. 



This shows that 



(3.4) 



Af-l 

E 



where 



RN{r, u) 



^ (7v-i)!r(i.) 

(7v-i)!r(^.) Jo 



{1 + v) 



r u / —Q Q 



N-1 



-dv 



[r "m" — v) 
I 1 + e~2~w j 



dv. 



Notice that 



Therefore 
(3.5) 



1 + e 2 w 



u+N 



< 



1 



u-\-N na 



sm 



\RN{r,u)\ < 



r(N + i') 



Nir{iy) sin 



h'-\-N TTCt 



Substituting p.4p into (|3.3p . we find that 



(-f)^r(j. + j) 



2 2 TT 2 



J=0 



+ / K n-2 (u)RM(r, u)u 2 (iM 



The assertion (|3.f p follows from the formula (see e.g. [^) 



Kr^{u)u'^+'^Uu = 2t+"J'-ir 



(3.6) 

and the bound p.Sp . 

An immediate corollary of Proposition 13.41 is 



□ 
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Corollary 3.5. When a G (0,2) and v > 0, the asymptotic expansion of qa^v.nir) 
as r —!■ oo is given by 

(3.7) 



TT— —J j! r(j/) 

For the special cases where v — \ ov n — \, (|3.7p was obtained in [551 [Ml Ell HH 
[13] by different methods. Notice that the large r-asymptotic behavior of qa,u,n{f') 
is very different for the case of a e (0, 2) and a — 2. qa,v,n{'r') decays polynomially 
when a £ (0, 2) and exponentially when a — 2. Also notice that the large-r leading 
term of qa,u.7i{T) is of order r~"~", which is the same as the large-r leading term 
of Sa,„(r) iMI). 

Now we turn to the asymptotic behavior of qa.y,n{r) when r ^ 0. For this 
purpose, we need to derive an integral representation of qa.y,n{'i') which extends the 
results of [11[S51[II1[I11[I3]- We begin with the following lemmas. 

Lemma 3.6. Let z he a complex number such that \arg{z)\ < n, then 

1 1 /-c+ioo 

where < c < v. 

Proof. Using the fact that when u 3> 1 (see e.g. [3]), 

(3.9) \T{c + iu)\ = V2^\u\^-h-^ (^l + o(^-^^ 
we find that if u ^ 1, 

(3.10) r(c + iu)r(i^-c- w)z-^-^" < C|zr'=el"ll''''S"l|u|''-ie-''l"l. 

Here and the foUowings, C or Ci, C2 represent constants whose values can be differ- 
ent in different lines. p.lOp shows that the right hand side of p.Sp indeed defines 
an analytic function of z when |arg(z)| < tt. Now using 



we have 

1 1 



(3.11) 



{1 + zr T{v)j^ 



— / r{w)z-dw, 

^" ' J c—ioo 



1 / u~l-y I ^ 



T{v) Jq [_ 2TTi 

For u > 1, ([3^ implies that 



y^^'e-y — / T{w)y'-z--dw dy. 



Therefore we can interchange the order of integrations in (|3.1ip when |arg(z)| < 
7r/2 to obtain p.Sp . The general case where |arg(z)| < tt follows by analytic 
continuation. □ 
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Lemma 3.7. For a E (0, 2) and u > 0, we have 

1 11 rc+ioo 

Im- : — ^ ——— Y{w)r{iy -w)u~"'^ sin— -dw, 



where < c < v. 

Proof. From Lemma 13.61 we have 

Im -Tj =Im <{ T-^-=rr^ I T{w)r{i^ - w)e-'^u-~"'^dw 




, , , r 



c+'ioo 

r 

2TTi r(i/) J^_i^ 



{w)Tiiy - w) cos (^) u-"'"dw 
{w)r{iy ~ w) sin (^) ^""'"du; 



Therefore, it suffices to show that 



27rz r(z/) 



C — lOO 



and 

Notice that the functions z i— > sinz, z i— > cos z, z i— *■ u^, u S M and z l—^ T{z) are aU 
real on the real axis. Reflection principle implies that they satisfy /(z) — f{z). On 
the other hand, 

V^lh r^^T{w)nv-w)coJ^)u~--dw 

By taking the complex conjugate of p.l4p . it is easy to see that p.l4p is real. 
Therefore p.l2[) holds. (|3.13p follows analogously. 

□ 

Now we can prove the following. 
Proposition 3.8. The function qa.v,n{i') has the following integral representation: 

j.-n 2 l-c+iw 

(3-15) qccuA'r) r,^, /Q,,.,n(w;r)dw;, 

where < c < minji^, — }. Here 

(3.16) f^^,Aw;r)=T{w)ri^-w)[-j -fj^- 

Proof. Using (|2.10p and Lemma [5T71 we have 
(3.17) 

qa,u.n{r) = „^-2 „+2 / (™)ut <^ T^T:rr^ / r{w)riiy - w) sin ZL^u" 

2 — TTT- Jo " 1 27r« r(:/) 2 
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For < u < 1 and ?/ > 1, (j2J2|) and (HH) imply that 
(3.18) 

2 2 

For M > 1 and y > 1, ([^1^ and dSH) give 
(3.19) 



Since c < I /a, p.lSp and p.l9p show that we can interchange the order of inte- 
grations in (|3.17p and use l\3.6\i to obtain 

^--^^^ ^ r(Hr(.-.)sin^ (-) r (^^j r 

()3.15p follows by using the identity 

r(z)r(i-^) 



SlUTTZ 

□ 

By making a change of variable, it is easy to see that (|3.15p coincides with the 
results of [IHl ISH [HI [HI US] in the special cases where v — 1 or n = 1. However, 
instead of proving that the right hand side of (|3.15p satisfies the equation (|2.2p . we 
choose to give a direct derivation of (|3.15p here. 

To find the asymptotic expansion of qa^v,n{T) as r 0, we can apply residue 
calculus to (|3.15|) . More precisely, we shift the contour of integration to the right 
and pick up the residues of fa,v,n{w]f)- Since the function T{z) has simple poles 
at the point z = —k, k = 0,1, 2, ... with residue 

(3.20) Res,=_fcr(z) = 

the poles of fa,v,n{w; r) that are situated to the right of the line Re u; = c are at 
the points 

n + 2? 

w = -, j = 0,1,2,... and w = v + l, 1 = 0,1,2,.... 

a 

These are simple poles except if n + 2j = a{v + I) for some {j, Z) G N x N. Here 
N = N U {0} is the set of nonnegative integers. Define the set A„ by 



A, 



v) e (0, 2) X M+ I n + 2j = a{v + /) for some (j, e N x n| 



so that {a,v) g A„ if and only if fa.u,n{w;r) has double poles on the right half 
plane Re w > 0. The detail characterization of the set A„ is deferred to the end of 
this section. In the following, we give the asymptotic expansion of qa,v,n{T) when 
r^O. 

Proposition 3.9. When a G (0, 2) and v > 0, the following formula is valid: 

N{r), 
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where 
(3.21) 



for any j^N 



p / n-a(?y+0 



2 ^ ^_^T{^ fai^-n-2j_\rr_Y^ 



0<j< 



/or any /^N 

> -i-V x ; \ { -^b(v + + 



/or some j^N 
0<1<N-1 

-H 2 +v+^°g(2) jy 

and 

(3.22) i?a,.,.;iv(r) = O "+"(A^-|)) 



0. 



A^ote that the third summation in p.2ip is zero if and only if {a,u) ^ A„. 

Proof. We use the integral representation of qa.^,n{r) given by (|3.15p . Given a 
positive integer N, we observe that there is at most one point of the form w ~ "^^^ ; 
j £ N in between i/ + N — 1 and i' + N since a < 2. If there is no such point, let 
Qn = 1^ + N - ^. If WAT = , jAT e N, lies between u + N - 1 and ly + N - ^, 

let Qat = 1/ + iV - i. If WAT = "+^^" , Jn G N, lies between v + N - ^ and + iV, 
let Qat = iy + iV - |. Notice that 

(3.23) mm{\QN-{iy + N-l)\,\QN-{iy + N)\}>^, 
and if there exists wn lying between v + N — I and u + N , then 

(3.24) minjlQAT-u'Arl} > i 

Consider the contour Cn,m which consists of the boundary of the rectangle 

{«) I c < Rew < Qni |Imu>| < M}. The poles of fa,v,n{w] r) inside Cn,m are the 



points w = V + I, I ^ 0,1,. . . ,N - 1 and w = ^^^±21^ j = 0, 1, . 
When c <Rew < Qn and Imtu = ±M, ([SJ]) shows that for A/ » 1, 



a(i/+Af-i)-r 



|/a,,.,n(w;r)| < CAf*-"'=+''-ie-^*^2-"'=max{l,r}"(''+^). 
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Therefore, the integral of fa,v,n{w]r) over the two horizontal segments of the rec- 
tangle Cn.m goes to zero as M oo. This gives 



z inside C 



If (a, J^) ^ A„, all the poles of fa,u.n{w; r) inside Cn,m are simple and (j3.20p implies 
that 



2 inside CiN 



N-l 



(3.25) 



E 



-p / n— Q(f/+0 
V 2 



(-1) 



■T{u + l)- 



y+l) 



E Res^._,>+2j/„,,^^„(u;:r) 
2/ 



+ 



E 

J=0 



This gives p.2ip when {a^v) ^ A„. On the other hand, if {a^v) e A„, then there 
are (j, e N x N such that n + 2j = Q;(i^ + Z). In other words, / = - v for 
some j e N and j — "('^+')^" for some I e N. We have to omit such terms from 
the first and second summations in (|3.25p and add the residues of fa,u,n(w',r) at 
w = V + I = which gives 

(3.26) 



E 



(-i)'+ °"^+2"'" r(z/ + /) 



+')-" ^ I 



2 



for some j/'GN 
0</<iV-l 



2 ^ 



Q!(t^ + ^) — n 



1 



This gives the third term in ()3.2ip when ioi^v) e A„. To obtain ()3.26|) . we have 
used the formula 



r(-/c + z) 



fc! 



+ V(fc + 1) ^ + 0{z) as z 0, 



fc e N. 



Finally we consider the remainder term. For y large enough, (|3.9p gives 



(3.27) 
Therefore 



fa.jy,n{QN + iV\r) 



1 
2^ 



Qjv+ioo 



Qjv — ioo 



which gives (IS22)- 
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□ 

By setting n — 1 or ly — 1, one can check that the resuhs of [TTl [T^] and [52] 
were recovered respectively. From (I3.2ip . we see that qa.v,n{r) ^ oo as r ^ if 
and only if av < n. More precisely, we have 

Corollary 3.10. Let a E (0, 2) and v > 0. As r ^ 0, the leading term of qct^^.n{r) 
is given by 

/. If ai^ < n, then qa.v,n{f) • i ^ 

II. If av > n, then qa,v,n{i') 



1 ri 



2«-ia7rtr(i/) r(f)" 

///. Ifav = n, then g^.^.^W - |- log (0 ' + V (f ) + ^(1) - ^ - ^(1)) 



One notice that putting a = 2 in the results of Corollary 13.101 one obtains the 
results of Corollary 13.31 We also remark that the large-r asymptotic expansion of 
qa,u,n{f') when a E (0, 2) (|3.ip can also be obtained from (|3.15p by considering the 
poles on the left of the line Ke w — c. 

Now we analyze in more detail the conditions for which (a, i^) G A„. Observe 
that it is necessary that both a and are rational numbers or both are irrational 
numbers. If a and v are both irrational, then 

n + 2ji = a{iy + h) and n + 2j2 ^ a{iy + for (ji, ^i), (ja, ^2) e N x N 

if and only if 

2(ji-i2)-aGi-Z2), 
if and only if ji = j2 and h ~ h- In other words, if a and v are both irrational 
numbers, there is at most one pair of nonnegative integers (j, I) such that n + 2j = 
a{i' + /). For the case where both a and f are rational numbers, we write a = a/b 
and ly — e/ f, where a, b,e, f E N, gcd (a, b) = 1 and gcd (e, /) = 1. Then 

(3.28) n + 2j = a{iy + l)=^n + 2j = ^^^-j^, (j, e N x N. 

Since gcd (a, 6) = 1, gcd(e,/) = 1 and the left hand side is an integer, it is 
necessary that / divides a and b divides e + If. Let a = mf for some m E N. Since 
gcd (a, 5) = 1, this implies that gcd (&, /) = 1. Therefore, it is always possible to 
find positive integers k and / such that 

(3.29) e + lf = kb. 
We have then 

a{e + lf)_ 

However, there does not necessary exist a nonnegative integer j such that n + 2j = 
km. We need to discuss the parity of n. If n is odd, then n + 2j is odd for any 
j E N. Therefore for n + 2j = km, it is necessary that m and k are both odd. 
TO is odd implies that both a and / are even or both are odd. In the case where 
a and / are both even, b and e are both odd. Then for any I E N, e + If is 
always odd. Therefore if {k,l) G N x N is a solution to p.29p . k is necessary odd. 
This implies that there are always solutions {j,l) g N x N to p.28p . Moreover, if 
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{jo,lo) is the smallest nonncgativc solution, then all the other solutions are given 
by (jqjq) — Uo + ag/2, Iq + bq), g G N. In the case where a and / are both odd, if 
{k, I) is a solution to ()3.29|) . then so is {k + f,l + b), but one of k and k + f must 
be odd. In this situation, we also find that there always exist solutions to (|3.28l) . 
Moreover, if (jo, Iq) is the smallest nonnegative solution, then all the other solutions 
are given by (jg, Iq) = {jo + aq, Iq + 2bq),q e N. 

For the case where n is even, n + 2j is even. In this case, m can be either odd 
or even. If m is odd, k has to be even. The same reasoning above shows that 
a and / has to be both odd and there is always solutions to (j3.28p . Moreover, 
if (ioi^o) is the smallest nonnegative solution, then all other solutions are given 
by (jqilq) = (jo + a?, ^0 + '2bq),q G N. If m is even, then a is even and there 
is no restriction on A:. In this case, there is always solutions to (|3.28p and all 
the solutions can be expressed by the smallest nonnegative solution (joj^o) via 
Uq, Iq) = Oo + aq/2, lo + bq), qeN. 

We summarize the results as follows: 

Proposition 3.11. Let a E (0, 2) and i' > 0. 

I. If a and v are both irrational numbers such that {a,v) £ A„, then there is a 

unique solution [j, I) to the equation n + 2j — a{i' + I), [j, Z) G N x N. 

//. If a and v are rational numbers with a — a/b, v — e/f, a,b,e,f £ N and 

gcd{a,b) = gcd{e,f) — 1, then (a, i^) G A„ if and only if f divides a and one of the 

following holds: 

(i) n is odd and a/ f is odd. 

(a) n is even and if a/f is odd, then a and f are both odd. 

Moreover, if{jo,lo) is the smallest nonnegative solution ofn+2j = a{h'+l), {j,l) G 
N X N, then if a is odd, all solutions are given by {jq,lq) = [jo + aq, Iq + 2bq) , q G N. 
If a is even, all solutions are given by (jq, Iq) = (jo + CLq/2, Iq + bq), g G N. 



4. Representation of qa,i,,,iir) in terms of entire functions 

In this section, we are going to analyze under what conditions the asymptotic 
expansion of qa.v,n{'r) given in Proposition 13.91 will give rise to representation of 
qa.u,n{T') in terms of absolutely convergent power series. First we have 

Proposition 4.1. Let a G (0,2) and v >Q. Then 

qa.u,n{r) = lim Sa.u,n-N{r) 

N^OG 

uniformly for r in any compact subsets of (0, oo). Here Sa,i^,n;N{'r') is given by 
(13^11) . 

Proof. Following from the proof of Proposition l3.9[ the remainder term Ra,v,n;N {f) = 
(r) is given by 

1 rQN+ioa ^— n ^ 

Ra,u,n-N{r) = ^rq-r— fcc.u,n{w\'r)dw = — 

TT 2 T{v) 2m Jq^-ioo tt 2 r{v) 2tt 



fa,u,n{QN + w;r)dy. 
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where fa,u,n{w]r) is given by (|3.16|) . We rewrite it as 

(r/2)"('5iv+»a)r(Q^ + iy) 



fa,i',n{QN +iy\'r) 



(4.1) 



r(l - J. + + iy)T |^ 2+a(Q»+»y)-n ^ p 



Sin TT (Qjv + ~ v) sin 



We want to find an upper bound for \ fa,i^.n{QN + iy',r)\ when Qn is large. Since 
fa,u,niQN + iy;r) = fa.,u,n{QN - iy;r), it suffices to consider y > 0. By Stirling's 
formula (see e.g. [5]), 

(4.2) r(z) %/27rz^"3e-^ as \z\^oo, |argz| < tt. 

Therefore, if < v < u, 

Ciw"-5e-"-T" < \T{u + iv)\ < C2(%/2w)"-5e-"; 
whereas if w > w > 0, 

Ciu"-5e-«-t" < |r(u + w)| < C2(%/2w)"-^e-"-^^ 
These imply that if Qat 3> 1 and < y < Qn, 



r(QAr + iy) 



r(l - + Qat + iy)r (H±£i 
and if Qn < y, 



(QN+iy)- 



(QN+iy) 

2 



< C 



(V2(2e/a)"e("+i)^)' 



r((9A' + iy) 



2+a(QN+iy) — J 
2 



r(i - i/ + Qat + iy)r ( 

On the other hand, if u > 0, we have 



^(QN+iy) \^ 



< c- 



(V2e") 



Q 



j+aQjv + i — I'— -S- 



-y 



N 



1 




2 


sin(u + iv) 







< 



e" - e 



Therefore, if 1 <C < y, 

\fa.v,n{QN + iy\r)\ < C- 
( 

and it follows that 



(\/2e")'5" 



1 371 / T* \ 



,Qjv+qQ]v + ^-i'--| 



-y^" 2e 4 



(9 



(Q A + W,r)dy 



\y\>QN 



<C 



(%/2e") 



(y2e""i)'3" /r\°'Qf 



~ ^aQjv + i-i^-f V2/ 

which tends to zero when N — > oo, uniformly for r in any compact subset of [0, oo) 
For y < Qn, p.23p implies that 



1 



sin7r((5jv -\- iy - v) 



< 



sin7r((57v — i^) 



< 



1 



18 



S. C. LIM AND L. P. TEO 



and 



Consequently, 



i 



l^/l<Qi^ 



1 

T(a(QN+iy)-n) 



fa.p,n{QN + iy;r)dy 



< 



< C 



1 

7r(aQiv~K) 



< 



which tends to zero when N oo, uniformly for r in any compact subset of 
[0,oo). □ 

Now we want to investigate the analyticity of the series representation oiqa,iy.n{'r')- 
First we have 

Proposition 4.2. Let a = a/b G (0,2) and v = e// > fee rational numbers with 
gcd(a,b) = gcd(e,f) — 1. If one of the conditions 
(i) a is not divisible by f , 

(a) n is an odd integer and a/f is an even integer, 

(Hi) n is an even integer, a and f are both even integers and a/f is an odd integer, 
holds, then qa.v,n{'r) can be represented by: 



(4.3) 



Q_a.,y.n (^) 



2"7r-r(!y) 12 
where Ai{z) and ^2(2) are entire functions given by 



Ar;.,,„(z)=5:^l^r(. + 0-^ 



r 



(71 — a 



(4.4) 



i=0 



r 



^{y+i) \ 2" 



a 



22j 



Proof. If the conditions (i) or (ii) holds. Proposition 13.111 implies that {a,v) ^ A„. 
Using Propositions 13.91 and 14.11 we only need to show that the right hand sides of 
the equations ()4.4|1 define entire functions. Let 



(-1)' ^ 



l — Q.{v-\-l) 



(-1) 



T{v + l) 



II 



r 



2 



2,i 



(-ipr(^) p ^a:. - n - 2j^ _ (-1)^-1 r (I^) 
j! r(^) 



Given Z e N, there exists a unique ji such that 



a(i^ + Z) — n 



Ji 



< 



Since 



2 jl + 0, 

Q!(i^ + /) — n 



31 



a(e + fl)-bf{n + 2ji) 
26/ 



> 



26/- 



MULTIVARIATE GENERALIZED LINNIK'S PROBABILITY DENSITIES 



19 



Similarly, given j e N, there exists a unique Ij such that 



n + 2j — av 



< 



n+2j — ai/ 



Ij ^ implies that 



n + 2j — av 



-I. 



bf{n + 2j) — ae — alj 



af 



> 



Therefore, 



sm 



7r(a(z^+/) — n) 



Lvr I — 



sm vr I — ^ — ji 



< 



sm 



2bf 



and 



1 



sm 



n(7i-\-2j — au) 



< 



sm ^ 

af 



3m ^( i I,) 

Together with Stirling's formula (14. 2p , we find that for I and j large enough. 



i+i| 



1^1. 



< 



1°" sin 



and 



2fc/ 



< 



2jl 



7^ sm — f 



Therefore, the two infinite series on the right hand sides of the equations ()4.4|) 
converge absolutely and uniformly on any compact subsets of C. □ 



Proposition 4.3. Let a = ajh S (0,2) and v = ej j > be rational numbers with 
a, 6, e, / S N and gcd{a, b) — gcd{e, /) = 1. // one of the conditions 

(i) n is an odd integer and a/ f is an odd integer, 

(ii) n is an even integer and a/ f is an even integer, 

(Hi) n is an even integer, a/ f is an odd integer and a and f are both odd integers, 
holds, let (joj^o) be the smallest nonnegative integers satisfying n + 2j = a{v + 
0,(j,0 e N X N. Then 

Qa,i^,n(^) can he represented by: 



Qa.^u,n (^) 



1 



(4.5) 



2"7r-r(i^) 1 2 

+ 22l^A3;a,.,„(r'^) + 22j^A4.a,.,„(r'^)l0g- ^, 
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where if a is odd, Ai{z), A2{z), ^3(2) and Ai{z) are entire functions given 
(4.6) 



leN: l^lo+2bq,qeN \ 2 

" ^ ('o + 2g6)!(jo + ag)!r(f+jo + ag) [ 
- i/'(^o + 2bq + 1) - (I + io + ag) - ^V' (io + + 1) 



22a9 ' 



A, (z) 2f-l)'o+^o V (~1)^ r(. + ^o + 25g) 

A4;.,..„(-)- 2( 1) Z.(;^ + 2g6)!(jo + ag)!r(t+io + a9)2^''^' 

whereas if a is even, 
(4.7) 



(-l)'w.. , n^l^ 



n— a(i>'+0 



^i-Qi/n(2)= . r(;^ + 

i,Q,i.,nV ; ^ n ^ ' p/a(iy+0\ 2"' 

;eN; l^la+bq,qeN V 2 ^ 

,^ 2(-l)'o+Jo- (_l)9(''+f) r(i. + ;o + fog) , 

^^^"''^^"^^^ ^ ^ g (/o + .5)!(.o + f)!r(t + ,„ + f)\^^(^ + + 

A4;....„(^)- 2( 1) Z.(;^+^^),(^.^^^),r(n+,„ + ^)2a,- 

Proof. The fact that 9Q.i/,n(?') can be represented in the form (|4.5p foUows from 
Propositions 13. 9i 13.111 and 14.11 The proof of the absolute convergence of the series 
Ai:a.iy,niz) and A2-^a.u.n{z) , ^3:Q.i/,n (-z) and ,y „(z) foUows the same as the proof 
of Proposition [121 and the formula (see e.g. [5]): 



Tpix) — logx — - — y O [ ^ ] asx^oo. 
2x \^ J 

□ 

Propositions 14.21 and 14.31 show that when a and v are both rational numbers, 
qa.u,n(r) can be represented in terms of entire functions. In the case where a is a 
rational number and v is an irrational number, we have 
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Proposition 4.4. Let a ~ a/b E (0,2) be a rational number with a,b E N and 
(a, 5) = 1. If v > is an irrational number, then qa,u,n{f') can be written as the 
sum oj two absolutely convergent series as in Proposition] 



Proof. Notice that the conditions a is rational and v is irrational imply that {a,v) ^ 
A„. Given ^ > an irrational number, let be the unique positive integer such 
that 



a 



This implies that for all g G Z 



(4.8) 



> min < (5 



< 
1 

2a 



1 

2a' 



Given I g N, there exists a unique ji such that 

a{v + I) — n 



31 



1 

< -. 

- 2 



implies that 

a{v + I) — n 



31 



b{n + 2ji)-al 



> 



Similarly, given j E N, there exists a unique Ij such that 



2j 



< 



1 



(|4.8p implies that 



n + 2j 



I, 



The rest of the proof is the same as Proposition 14.2 



' + 2j) — alj 



> 



□ 



When is a rational number but a is irrational, the situation is more compli- 
cated. We have to introduce the concept of Liouvillc numbers as in [551 IHH [HI 
[T2l [T3] . Recall that an irrational number /3 is called a Liouvillc number if for all 
TO = 2, 3, . . ., there exists a rational number p/q,p, q E Z, such that 



q 



< 



By Liouvillc theorem (see e.g. |63|). all Liouvillc numbers are transcendental and 
the set of Liouvillc numbers, denoted by L, has Lebesgue measure zero. 

Proposition 4.5. Let v = e// > &e a rational number with e, f E N, (e, /) = 1. 
If a is an irrational number but not a Liouvillc number, then qa.v,n{^) can be written 
as the sum of two absolutely convergent series as in Proposition\4-.2\ 



Proof. Notice that the conditions on a and ly imply that (a, t^) ^ A„. Since a is 
not a Liouvillc number, there exists an integer such that for all m > TOq, and 
for aU p e N, g e N, 
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Fix I G N. There exists a unique ji e Z, such that 



a{iy + I) — n 



< 



By (|4?9)) . we have 

a{i' + I) — n 



31 



{v + l) 



a 



f{n + 2ji) 



e + fl 



> 



This imphes that 



(4.10) 



1 



T{a{u-\-l) — n) 



( a(y^l)-7i ■ \ 

smni^^-^ Ji) 



< 



2/™°(z/ + Z)""-i 



On the other hand, given j E N, let Ij be the unique integer such that 

1 1 (n + 2j) , , , 1 1 
a 2 a 2 a 

This imphes that 



1 



> 



v + lj n + 2j + 1 



Consequently, 

(n + 2j — av) 



a 



f{n + 2j) 



> 



e + fl, 
1 



> 



and 
(4.11) 



1 



sm 



7T{n+2j — au) 



< 



/™°(n + 2j + l)™»-i ■ 
/™°(n + 2j + l)™°"^ 



TTa" 



If iTia is larger than a, the ratio jAi^j+i |/|v4i_; | or |yl.2j+i |/|v42j | may become un- 
bounded as Z ^ oo or j oo. Therefore, instead of considering the ratio, we 
consider the roots \Ai^i\t and |A2j|j. Stirling's formula (|4.2p . (|4.10p and (|4.1ip 
imply that 



\All\T <Cl 







and 



-f — ^ 



0. 



This shows that both the series Ai.a,u.n{z) and A2-a,u.n{z) are absolutely and uni- 
formly convergent for all z G C. 

□ 



Recall that an entire function E{z) = J2T=o -^kZ^ is said to be of order p if and 



if 



k log k 

p^limsup- . 

fe^oo logl/l-Efcl 



It is easy to deduce from the proof of Proposition 14.51 that when a, v are both 
rational numbers, or if one of a or is a rational number, and the other is a non- 
Liouville irrational number, then the entire functions Ai-^a,v,n{z) and A2-a,v,n{z) 
have order 1/a and 1/2 respectively. In the case where a = a/b and v = e// satisfy 
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the conditions (i), (ii) or (iii) of Proposition [473l the entire functions A^-a^v^niz) and 
A4-a,iy,niz) have order 1/a. 

Although Propositions 14.21 14.31 and 14.41 are proved for a G (0,2), one can show 
that by naively extending the results to a = 2, one obtains exactly the results for 
a = 2 given in Proposition [3?2l 

At this point, one may wonder whether there exist any values of a and v where 
both the series ^i;Q,;y.n(z) and A2;a.v,n{z) are divergent. As was proved in [62], 
when V = \, there is a dense subset of a S (0, 2) where both the series Ai.^a,v,n{z) 
and A2-a,u,niz) are divergent. Now we extend the result to general > which can 
also be regarded as the generalization of a result in 12J for n = 1. 

Proposition 4.6. Given a rational number v = e// > 0, there exists a dense 
subset of a €z (0,2) where both the series Ai-^a,v.n{z) and A2-a.v,n{z) are divergent. 
This implies that there is a dense subset oj {a^v) e (0,2) x (0,oo) where both the 
series Ai.a.v,n{z) and A2-a,iy.n{z) are divergent. 

Proof. The proof follows closely the ideas of [62l [12]. Given ly = e/f, define a 
sequence of positive integers {cfe}^! by ai = 1 + 2/ and ak+i = (1 + 2/)^'^'=. Let 
A be the set of sequences {^m}m=i satisfying the conditions: 

(i) Sm is equal to or 1. 

(ii) Sm = 1 impHes that m G {(^k}kLi- 

(iii) 6m — 1 for infinitely many m. 

Define 17 to be the set of real numbers of the form 

C30 ^ 

and let A be the set of rational numbers of the form 
t 



(1 + 2/)-^ (1 + 2/)- ' s,ieN,a„, e{0,l,2,...,2/}. 



m=0 

Finally define E to be the set 

E ^ ia e {0,2) : av ^ — = x + y, x e A.y e fl 

Clearly i? is a dense subset of (0,2). Given a ^ E, we claim that for any r > 0, 
there exists a subsequence of positive integers Ik and a subsequence of integers jk 
such that 

|^i,ifcz''°| ~*°°> oo and \A2^j^z-'^\ ~*°°> oo 

for any z ^ 0. This will imply that the series Ai-^a.v,n{z) and A2-a.v,n{z) are 
divergent. Given a € E, av can be written uniquely as 

t CO c 

av = h+y , , + y , , , /leN, a,„e{0,l,...,2/}. 

m=l ' ' m=t+l ^ •' ' > 

Define the sequence {r]k}'kLi so that {m > t : Sm = 1} = {ilk}kLi- Notice that 
rji > t + I and 

Vk+i > (2/ + 1)2"- >9r]k. 

For any fc e N, we have 

( t k ^ \ oo ^ ^ 

^ < 



' r ^ £ (2/ + 1)™ ^ S (2/ + 1 (2/ + D- ^ (2/ + 1)--- 



24 



S. C. LIM AND L. P. TEO 



This implies that 



aj.(2/ + 1)"'' - h{2f + + ^ a,„(2/ + l)*"'" + ^ (2/ + l)'""""" 



< 



< 



(2/ + lyik+i-vk-i (2/ + l)ik+i/2 ■ 
Notice that since (2/ + 1) is odd, the integer 



Pk = h{2f + 1)'"= + J2 am(2/ + 1)"'="™ + E (2/ + 1) 

m— 1 m— 1 

is odd if and only if 

t 

h + E flm + fc is odd. 



Vk-Vn 



m — 1 



If ft, + X]L=i — is even, p2fe — f^ is even. In this case, we define Ik and jk by 



e [(2/ + 1)"^'= - 1] 
/ 



Jk 



P2k - n 



If ft + J2m=i ^rn — n \^ odd, P2k-i — n IS cvcn, and we define Ik and by 



:[(2/ + l)'^2.-i _l] 
/ 



P2fc-i - n 



In the first case where ft + X]L=i — n is even, we have 



+ lk) - n 



■ Jk 



au{2f + 1)"^'= -n-p2k + n av{2f + 1)"='= - p^k 



< 



2(2/+ l)''2fc+i/2 



and 



n + 2jk - av 



a 



a[v + Ik) — n 



■ Jk 



These imply that 



1 



1 



+ifc) — w 



Jk 



< 



> 



>;(2/ + l)'te'=+i/2' 



2(2/ + 



and 



7r(Ti+2jfc-at/) 



1 



smTT 



+ifc)-n-2jfc 



> 



a 



(2/ + l)')2fc + l/2 



Together with Stirling's formula (|4.2p . we find that for fc ^ 1, 

, C(2f + l)'"=+i/2 , C(2f + l)''2fe+i/2 
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Notice that h < v{2f + 1)''^'= and jk<{h + t + k){2f + l)"^^" . This impUes that for 
k large enough 



(^(2/ + 1)3...) 



(2/ + !)"-'°+^/^ _ :ii|ti_4.,,,(2/+i)... 

-((2/ + l)4')2.)-(2/+l)''='= ^^-^^'^ 



Now 



^?H|ti - 4iyr,2k{2f + 1)"- >^(2/ + l)'"^" - 4i.r;2fe(2/ + l)"^-" 

=(2/ + 1)"- (^I^^il!!^ - 4z.r?2.) ^ oo. 

This shows that \Ai^i^z^'' | ^ oo as fc ^ oo. The proof for \A2j^^z^'' | ^ oo as A: ^ oo 
is similar. For the case where h + a™ — n is odd, the proof is the same. This 

concludes the assertion of the proposition. □ 
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